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Abstract 

We investigate the low temperature asymptotics and the finite size spectrum of a class 
of Temperley-Lieb models. As reference system we use the spin-1/2 Heisenberg chain with 
anisotropy parameter A and twisted boundary conditions. Special emphasis is placed on the 
study of logarithmic corrections appearing in the case of A = 1/2 in the bulk susceptibility 
data and in the low-energy spectrum yielding the conformal dimensions. For the sZ(2|l) 
invariant 3-state representation of the Temperley-Lieb algebra with A = 1/2 we give the 
complete set of scaling dimensions which show huge degeneracies. 
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1 Introduction 



Many different, and at first sight unrelated physical phenomena like the thermodynamics of 
quantum spin chains, critical properties of two-dimensional classical vertex models and of su- 
persymmetric network models are related on fundamental mathematical grounds. Often, such 
relations are established by mappings of quantum systems to classical systems, where the related 
objects are operators like the Hamiltonian and the transfer matrix 'living' in the same Hilbert 
space of physical states. Sometimes, the relations are less obvious as the related objects are de- 
fined in rather different spaces, but turn out to be equivalent representations of some underlying 
algebraic structure. The structure relevant to our investigations is the so-called Temperley-Lieb 
algebra [1-3]. 

The seminal model for correlated quantum many-body systems and the standard reference 
model of the Temperley-Lieb algebra is the well-known Heisenberg model. In one spatial dimen- 
sion the system with nearest-neighbour exchange is exactly solvable for the spin-1/2 case [4]. 
Despite the integrability a comprehensive theoretical treatment is still missing. A notorious 
problem is posed by the calculation of correlation functions even at zero temperature, not to 
mention the corresponding properties at finite temperature. However, the asymptotics of the 
correlations in particular in the ground-state and at low temperatures are quite well understood 
by a combination of Bethe ansatz [1] and conformal field theory [5,6]. At exactly zero tempera- 
ture the spin-1/2 Heisenberg chain shows algebraic decay in its correlation functions and thereby 
constitutes a quantum critical system. 

Recently, the special point A = 1/2 of the anisotropic version of the Heisenberg chain, the 
so-called spin-1/2 XX Z chain, has attracted strong interest as certain mathematical simpli- 
fications occur in the construction of the eigenstates [7-9]. Also its classical counterpart, the 
six-vertex model with suitable boundary conditions, allowed for new insight into many combi- 
natorial problems like for instance the alternating sign matrices and boxed plane partitions, see 
e.g. [10,11] and references therein. 

In the case A = 1/2 like for other special points, the so-called roots of unity of the XX Z 
chain, a crossover of critical exponents leads to logarithmic corrections to the dominant critical 
behaviour. A notable case is the isotropic model A = 1 with logarithmic corrections to the low- 
temperature asymptotics, in particular for the susceptibility [12]. For the case A = 1/2 strong 
logarithmic corrections to the boundary susceptibility were found in [13]. 

In this paper we report on a low-temperature analysis of the A = 1/2 bulk properties that 
reproduce and confirm some of the findings in [13,25]. The main motivation of this analysis 
is the applicability to the study of all low-lying excitations of the XX Z chain with arbitrary 
twisted boundary condition allowing to calculate the energy levels beyond the conformal field 
theoretical results. 

The feasibility of such calculations is interesting because of applications to statistical mechan- 
ical models of the spin quantum Hall transition [14-16]. The localization-delocalization exponent 
was calculated in [16] from the mapping onto the two-dimensional percolation problem. The un- 
derlying system is a supersymmetric realization of the Temperley-Lieb algebra. Another model 
for describing the spin quantum Hall transition was introduced in [15] and investigated exten- 
sively in [17], however, still leaving open important questions. One of those questions concerns 
the relation of the two models. 

For the network model [16] the complete classification of the excited states resp. the scaling 
dimensions was performed by a Coulomb gas analysis in [33]. A combinatorial analysis based on 
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a decomposition of the partition function for the Potts model on finite tori in terms of generalised 
characters can be found in [34]. 

Here we present a systematic representation theoretical analysis of the eigenvalues of the 
Temperley-Lieb Hamiltonian corresponding to the network model on finite chains with periodic 
and 7r-twisted boundary conditions. Our analysis is based on [35] where a Temperley-Lieb equiv- 
alence with particular emphasis on the boundary conditions was established for vertex-models. 
The precise mapping of the generic irreducible sub-representations of Temperley-Lieb Hamilto- 
nians with periodic boundary conditions to the Heisenberg chain with suitable twisted boundary 
conditions allows for the determination of critical properties. Also in [35], the multiplicities of the 
sub-representations were worked out. The gained knowledge of the complete spectrum allowed 
for the computation of the thermodynamics of the Hamiltonians. 

In Section 2 of this paper an approach to the thermodynamics of the Heisenberg chain on the 
basis of non-linear integral equations for just two auxiliary functions is analysed [21] (see also the 
related approaches by [22-24]). For A = 1/2 an analytical calculation of the leading universal 
terms to the free energy and the magnetic susceptibility as well as the next-leading logarithmic 
corrections is presented. In Section 3 we study the critical exponents of correlation functions 
at zero temperature. The leading conformal finite-size spectrum of the Heisenberg chain with 
twisted boundary conditions and the next-leading correction terms in the case A = 1/2 beyond 
the CFT results are derived. In Section 4 we address the super symmetric network model and 
calculate the complete set of scaling dimensions as well as their multiplicities and logarithmic 
corrections. 



2 Thermodynamics and logarithmic corrections 

We investigate the thermodynamical properties of the antiferromagnetic isotropic spin- 1/2 
Heisenberg chain 

L 

h = 2jJ2lSjS* +1 + sys y J+1 + AS?S? +1 ] (1) 

3=1 

with anisotropy parameter A = cos 7 where 7(A = 1/2) = ir/3. 

We note that the elementary excitations ('spinons') have quasi-linear energy-momentum 
dispersion with velocity 

sin 7 3\/3 

v = Jtt >a=i/2 -rr J - ( 2 ) 

7 2 

The free energy per lattice site of the system is given by the following set of non-linear integral 
equations [21] for auxiliary functions a, and 21 = 1 + a 

v/3 f°° — 

logo(x) = ~ coshx +4> + J [k(x-y)log%.(y)-k(x-y-iir + ie)logW.(y)]dy, (3) 

and the corresponding equation for a, and 21 = 1 + are related to (3) by 'complex conjugation', 
in particular by exchanging a <H> a, and by changing the sign of </) f> — <p though being real as 
given by the external magnetic field <fi = 2(f-~/) ~^A=i/2 \fih. The integration kernel is defined 
by the Fourier integral 



1 r sinh U - 2] ffc 1 x 

K X ) = 7T 7 \ cos(kx)dk -> A =l/2 7T ~2—r — ■ ( 4 ) 

2 ~ J-00 2coshffcsinh (- - l) p 2jT ' sinh - r 
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Figure 1: Specific heat c(T) data versus temperature T (J = 1) for the spin-1/2 XX Z chain 
with anisotropy A = 0, 0.5, 1 and zero field h = 0. 



In terms of the solution 21 and 21 to the integral equations the free energy is given by 



These equations are readily solved numerically for arbitrary fields and temperatures by utiliz- 
ing the difference type of the integral kernel and Fast Fourier Transform. In Figs. 1 and 2 the 
results for the specific heat and susceptibility are shown for zero external magnetic field. At low 
temperatures the specific heat shows linear asymptotics c(T) ~ with velocity of the ele- 
mentary excitations v and central charge c = 1. The zero temperature limit of the susceptibility 
is Xo = l/(2w(vr - 7)). 

For A = 1 there are logarithmic corrections to the susceptibility xCH a * l° w temperatures as 
noticed in [12,20] and to the specific heat c(T). For a detailed analysis of the isotropic point we 
refer the reader to [12,25,28]. We like to point out that already for moderately low temperatures 
the correction terms are noticeable for x(T), see Fig. 2. 

There are other cases of the anisotropy parameter A with logarithmic corrections at low 
temperatures. One of these cases is A = 1/2 as noticed and analysed in [13,25]. Already a 
glance to Fig. 2 suggests that these correction terms must be of higher, i.e. less leading, order. 
Here we like to present a completely analytical low-temperature analysis (for h « T) for the 
case A = 1/2 from the non- linear integral equations. 

With view to the asymptotical properties of the driving term on the RHS of (3) it is conve- 
nient to introduce the following scaling functions 




(5) 



a £ (x) = a(-x-£), 



£ := ln{v/3) 



(6) 
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Figure 2: Susceptibility x(T) data versus temperature T (J = 1) for the spin-1/2 XX Z chain 
with anisotropy A = 0, 0.5, 1 and zero field h = 0. In the inset the low temperature behaviour 
of x(T) for A = 0.5 and 1 is shown. 



approaching well-defined non-trivial limiting functions in the low-temperature limit which satisfy 
logo+(x) = -2e~ x + 0(l/f3 2 ) +<P + ^(x) 

/oo . _ 
fc(x-y)log2l + (y) - k(x - y - ivr + ie) log 2lJ (y) dy, (7) 
-£ L J 

where ip£(x) accounts for the contribution of the functions 21 and 21 on the negative real axis 

/oo . _ 
k(x + y + 22)log^(y)-k(x + y + 2£-i7r + ie)log^(y) dy. (8) 
-£ L J 

Using these functions the free energy at low temperatures takes the form 

^ roo 



f = eo 



TTVf3 2 



e x 1oe 



2l^2l£2l £ 2l £ (x) cix + O (1//3 4 ) 



(9) 



In order to analyse the low-temperature asymptotics in detail we perform the following manip- 
ulation 



[(loga+)'(log2t+) - (log a+) (log *+)'] dx + 'conj' 



/oo 
[(_ 2e - + + ^ £ (x))'(lo g 2l+) - (-2e- s + + ^ £ (x)) (log 21+)'] dx + 'conj', 
-£ 

/OO / 4 \ 

(2e- x + Y> £ (x)')(log2l+)cfa - ./.log (l + e^J + 'conj', (10) 



where in the second line we have inserted (7) and used the symmetry of the integration kernel 
(kij(x) = kji(-x) where indices 1 and 2 refer to functions o, 21 and a, 21 respectively). In the third 
line we have performed an integration by parts with an explicit evaluation of the contribution 
by the integral terminals using the limits 



log aj(oo) 
log 0+ (00) 
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log2lJ(oo) = log +e + l^) 
±<l>, lo g a£(oo) = log(l + e-t*) 



(11) 



Next, we show that the LHS of (10) can be evaluated explicitly as it is a definite integral of 
dilogarithmic type with known terminals 



LHS 



2 / (log aj)' log 2lJeZx - log aj log 21+ 

-£ 



+ 'conj , 



= 2/ log(l + e z )dz- -01og(l + e^) + 2 / log(l + e z )dz + -^log(l + e^), 



7T 



(12) 



where the second line is independent of (/> (to be checked by differentiation) and the particular 
choice 4> = 'directly' leads to the last line. 

From the general expression (5) of the free energy / in terms of the functions 21 and 21 the 
detailed functional dependence of 21 and 21 on the argument matters even for the calculation of 
the asymptotic behaviour of / at low temperature. The usefulness of (10) lies in the fact that 
the integral of interest is related to another integral where only the asymptotic behaviour of 21 
and 21 at large spectral parameter enters. 

Combining (10) and (12) we arrive at a formula particularly suited for studying the asymp- 
totic behaviour of (9) 

/oo 2 a roo _ _ 

e-*log[2l+2Lt]d* = — + -<f> 2 - 2 J [^(z)' log 21+ + ^ c (z) , logOj| dx. (13) 

The first term on the right hand side determines the leading 0(T 2 ) term of the free energy, 
see (9), the remaining integral is a higher order term which can be calculated by explicit use of 
special properties of the integral kernel (4), especially the asymptotic behaviour 

1 



k(x) ~ —^xe 



This allows to approximate (8) by 

/oo 
k(x + y + 2£) log 
-£ 



2l£(2/)2l £ (y) 



dy, 



(14) 



(15) 



leading to a symmetric double integral 



/oo _ 
Y> £ (x)'log2l+(x) +^ £ (x)'lo g 2Lt(x) 
-£ L 

/oo /*c 
-£ i- 



2l+(x)2lJ(x) A;'(x + y + 2£)log 



2£ 



7T- 



-2£ 



(/: 



e -x log 2t+(x)2lJ(x) 



dx 



(/ 



%&(y)%z(v)\ d v dx 

e^log [2l-(2/)I £ (y)l dyV (16) 
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Amazingly, the integrals in the final version are those we calculated already in leading order 



a: 



e~ x log[2t£2lJ]d 



x 



1 /vr 2 4 - 

h -d> 

3 r 



(17) 



Hence 



e~ x log[2l+2lJ](ia 



e - * log[2l+2lJ]dx 



-- / [^(x)'log2t+ + V £ (x)'log2lJ 



dx, 



e log 



2l+(x)2l+(x) 



dx 



(/: 



e^log 2l-(y)2l £ (y) dy .(18) 



The final result reads 



/ e~ x log [21+21 c] dx 

J-s, 



2 ^ \ 2 
^ 4 ^2 

— i- -<r 

3 3 V 



(19) 



implying logarithmic corrections in T for the free energy and the zero-field susceptibility 



/ - e 



27TV 
7T 



7T 4 

Y + 3 ( 



£ = ln(v/T), (h«T), 



6v 



~ e - — T 1 + (20/tt) z 1 + 



T 2 ln(u/T) 
12^2 



1 + (20/tt) 



X 



47Tf 



1 + 



2T 2 In(u/T) 
3^ 



(20) 
(21) 



This result is in complete agreement with the effective Hamiltonian analysis of [13,25] where also 
higher order terms were calculated. In this approach, the occurence of the logarithmic correction 
stems from the interference of the first-order contribution of the TT perturbation to the Gaussian 
Hamiltonian with the second-order contribution of the cosine of the bosonic fields. The main 
motivation of our analysis is not the systematic extension of the low temperature corrections. We 
intend to compute the entire low-lying excitation spectrum of the quantum transfer matrix and 
by the (IhL correspondence the spectrum of the Hamiltonian up to order 1/L and (lnL)/L 3 . 



3 Critical indices and low-lying eigenvalues of the Heisenberg 
chain 

A major area of interest in statistical physics is the study of critical properties of models at 
phase transitions and, more specifically, the determination of the scaling dimensions x, from 
which the critical exponents can be deduced. For the critical two-dimensional six-vertex model 
and the associated spin-1/2 XX Z quantum chain a wealth of information has been collected by 
combinations of finite size studies of the eigenvalue spectrum and scaling relations by conformal 
field theory. 
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The energy levels of the corresponding quantum model on a finite chain with periodic bound- 
ary conditions scale with the system size L like (see [6] and references therein) 



2tt 

E x — En = — vx + o I — 



L 



L 



En 



TTV 

6L 



(22) 



where En is the ground state energy, E x are the energies of the low-lying excited states and v 
is the sound velocity of the elementary energy-momentum excitations. Equivalent formulae can 
be set up for the spectrum of the row-to-row transfer matrix of the two-dimensional model. 

For the critical six-vertex model the scaling dimensions have been treated by analytical as well 
as numerical methods based on (22). The dimensions x are given by the Gaussian spectrum [29- 
31]. The analytic method for calculating finite-size corrections based on linear integral equations 
for density functions of Bethe ansatz rapidities was applied to many models. However, this 
method has the shortcoming that it fails for the higher spin-S" XX Z chains and similar problems. 
There is an alternative procedure based on non-linear integral equations (NLIE) for suitably 
chosen auxiliary functions which proved to be applicable also in the hitherto problematic cases 
of higher spin-S* chains, i.e. for systems supporting filled seas of string configurations. At the 
same time the NLIE method proved flexible enough [32] to give a comprehensive analytical study 
of the six-vertex model spectrum. 

The eigenvalues are obtained from non- linear integral equations [21] for auxiliary functions 
o, and 21 = 1 + a 

x f°° — 

log a(x) = Llntanh — + C + / [k(x — y) log2l(y) — k(x — y — m + ie) log2l(y)] dy, (23) 

2 V-00 

and the corresponding equation for o, and 21 = 1 + a related by complex conjugation. These 
equations are very similar to those in (3). The only difference is the change in the driving term 



vj3 x 

- — OLln tanh — . 

cosh x 2 



(24) 



Note that the asymptotics of both functions are the same exponentials with different prefactors: 
— 2v(3 exp(— x) and — 2Lexp(— x). The constant C in (23) is iir<p/(ir — 7) where (p is the twist 
angle of the generalized periodic boundary conditions. (For certain excited states, the constant 
C may take different values.) 

The energy corresponding to a solution to the NLIE (23) satisfies an integral expression 
very similar to (5). Therefore the analysis of the finite size Hamiltonian can be carried out very 
closely to the analysis of the low-temperature analysis of the previous section. 

The main results for the general XXZ chain were obtained already in [32] and may be 
summarized for the \/L contributions to the energy and momentum eigenvalues as 



JL (£ - ie „) = -JL 



L 

2^ J 



27T2 



00 _ _ /"OO _ 

e^log 2t+2Lt(x) dx+ / e^log 2l £ 2T £ (x) 

00 /»oo 

Q~ x \og 2l+2Lt(x) dx- e~ x log Sl-SLTOr) 

£ L - 1 J-S. L 



dx 
dx 



(25) 



These formulae have a rather universal algebraic form as they are valid for all eigenvalues of the 
Hamiltonian with suitably chosen contours! 
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In [32] the leading contribution to the integrals was calculated 



e^log 2l±2lJ(x) 
-£ L 



dx = -2tt 2 



X((f)±s((f) _ l_ 

2 24 



(26) 



where 



*(v0 = xs , m {<p) = 1 -^s 2 + 2(i -L_ ( m - 1) 



s(<p) 



(27) 



are the 'scaling dimension' and 'conformal spin', and S 1 and m are integers. From these results 
with S = m = and (22) also the (effective) central charge follows 



c(<p) = 1 - 



6(y/vr) 2 

1 — 7/-7T 



(28) 



For the special case of interest (A = 1/2) we are able to calculate the next leading con- 
tribution to the integrals on the right hand side of (25). We note that (18) holds with the 
result 



/: 



e x lo2 



(x) I dx = it 2 [1/12 - (x(<p) ± s(<p))] + 7r 2 £e- 2£ [(1/12 - x(ip)) 2 - s 2 (ip))] , 

(29) 



where now 



£ = hiL^£e- 2£ 



InL 
IF' 



(30) 



These relations imply that the next-leading corrections over the (1/L) conformal field theory 
results for the energy levels are 



E x - E C X FT = -27TV [(1/12 - *(v?)) 2 - s 2 M)] 



InL 



(31) 



where v is the velocity of the elementary excitations. Our results agree with those of [25] and 
(slightly) generalize them to arbitrary twisted boundary conditions. 



4 The s/(2|l)-invariant network model: Temper ley-Lieb repre- 
sentation with A = 1/2 

Here we investigate the critical properties of the statistical mechanical network model derived 
in [14, 16] for the spin quantum Hall transition which may occur in disordered superconductors 
with unbroken SU(2) spin-rotation symmetry but broken time-reversal symmetry [18]. Using 
supersymmetry, and averaging over the random SU(2) matrices a classical 3-state model on the 
square lattice was derived with local Boltzmann weights resp. the i?-matrix given as a superpo- 
sition of the identity and the projector onto the sZ(2|l)-singlet state of nearest-neighbours. 

Thus, the study of the localization-delocalization transition was mapped onto the two- 
dimensional percolation problem. Furthermore, the local Boltzmann weights derived in [16] are 
a realization of the Temperley-Lieb algebra [1-3] and hence the model is related to some (inho- 
mogeneous) six-vertex model. For the symmetric case, the network model in [16] is related to an 
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integrable homogeneous six-vertex model. In our representation theoretical treatment we give 
the identification of the type of boundary conditions to be used for the six-vertex model which 
is essential for the determination of the scaling dimensions in a finite-size analysis of eigenvalues 
of transfer matrices and/or Hamiltonians. 

The Temperley-Lieb algebra TL(A) is the unital associative algebra over C generated by 
ei, e2, . . . , ex,_i with relations 

ef = A e u 

d ej ei = ei, if i,j nearest neighbours (32) 
ei ej = ej ei, if i, j not nearest neighbours 

depending on the complex parameter A = 2A, where A is the anisotropy parameter of the 
anisotropic Heisenberg chain which is the standard (faithful) representation of the Temperley- 
Lieb algebra. The periodic closure of TL(X) has one additional generator ei and relations (32) 
with 1 and L beeing considered as neighbouring indices. For a full description of the periodic 
(and more generally twisted) closure we refer the reader to [35] and references therein. 

According to the realization of TL(X) as a diagram algebra we use the following graphical 
notation for the operators ef. 



(33) 



where the two curved links correspond to the bra- and ket-states of the local sZ(2|l) singlet. 

Obviously, the parameter A appearing in (32) has the counterpart of a closed loop which 
for the model in [16] is 1 (= 1 — 1 + 1 due to the alternating standard and non-standard anti- 
commutation relations for fermions on odd and even lattice sites). These constructions may be 
generalized to the case of a sl{n + l|n)-invariant network model, again of Temperley-Lieb type 
with A = 1/2, but now with a total number of states per site d = 2n + 1. 

The statistical mechanical model of [16] has a staggered structure with two types of local 
Boltzmann weights alternating from row to row, cf. Fig. 3. The local Boltzmann weights can be 
given in terms of so-called i?-matrices which are defined by 

Rj,j+i( u ) '■= sin(7 — u) + sin(u)ej (34) 

where 7 = ir/3 and u takes certain values Qa and @b alternatingly from row to row. For 
general parameters &a,b the model is not exactly solvable. However, for identical interaction 
parameters 6 := 6,4 = @b, the row-to-row transfer matrices belong to a commuting family of 
transfer matrices defined on a square lattice with the same width and height as in Fig. 3 but with 
horizontal and vertical lines to which spectral parameters are associated. On the vertical lines 
we place alternatingly and O. Now consider a single row with spectral parameter v attached 
to the horizontal line. The corresponding transfer matrix is denoted by T(v) and consists of an 
alternating product of R{v — 0) and R(v — O) matrices, see Fig. 4 a). All matrices T(v) with 
arbitrary complex argument v commute due to the Yang-Baxter relation satisfied by (34). For 
v = the lattice in Fig. 3 is recovered. This is most easily understood by observing that R(0) 
is proportional to the identity operator and T(0) factors as shown in Fig. 4 b). The product of 
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Figure 3: Depiction of the Chalker-Coddington model. Each vertical zig-zag line corresponds to 
a local quantum space. Solid and dashed lines correspond to dual representations of sZ(2|l). In 
the general, inhomogeneous case, there are independent local interactions R(@a) and R(@b) on 
even and odd rows. 



two transfer matrices T(0), cf. Fig. 4 b), corresponds to two successive rows in Fig. 3 times a 
translation by two lattice units. For a lattice of square shape the translation operators reduce 
to the identity and the partition function of Fig. 3 is identical to Tr [T(0) L ] . 

Instead of working out the critical properties directly for the model depicted in Fig. 3 (with 
= Qa = &b we use the integrability of the transfer matrix T{v). Due to the commutativity, 
the eigenstates of T(y) are independent of v. The largest eigenvalue of T(0) is realized by the 
eigenstate with the lowest eigenvalue of the corresponding Hamiltonian H. Note that also the 
product T2(v) := T(v)T(y+@) is a commuting set of matrices generated by the complex variable 
v. For v = and by use of 'unitarity' for the i?-matrix, this matrix can be shown to reduce to 
the translation operator by two lattice sites. The logarithmic derivative of T^v) is nothing but 



The critical properties of the system, e.g. correlation functions of local operators, are described 
by the scaling dimensions which are obtained from the finite-size data of the low-lying excitations 
of this Hamiltonian with periodic boundary conditions or with twist angle it. In the dense loop 
representation of the model, periodic boundary conditions render closed loops to evaluate to 
1 (= number of bosonic states per site - number of fermionic states per site), for 7r-twisted 
boundary condition such loops evaluate to the total number of states per site d. 

As shown in [35] the sub-representations of the above Temperley-Lieb model appear as sub- 
representations in the standard reference Heisenberg chain, however with a new twist angle 
<p appearing. The representations are constructed by attaching A;- many singlets (k < L/2) to a 
pseudo- vacuum on a chain of length N = L — 2k. The dimension of the space Q P N of pseudo- vacua 



L 




(35) 



i=i 
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a) T (v) := 



b) T (0) : 





































o e o e o e 



c) t 2 (0) = 




Figure 4: a) Graphical depiction of a row-to-row transfer matrix T(v) with spectral parameter v 
associated with the horizontal line and alternatingly and associated with the vertical lines. 
Interactions are formulated for each vertex with local Boltzmann weight given by an i?-matrix 
R(u) where u is taken as the difference of the two spectral parameters of the intersecting lines. 
This construction renders the transfer matrix integrable. Note the factorization taking place for 
v = at every second vertex as R(0) reduces to the identity, a) and b). 



of length N for periodic boundary conditions is given by (see appendix) 



N 




dim(SP N )(d)= [- + !——) + \--^-—\ . ( 36 ) 



which specializes for a 3-dimensional local Hilbert space to dim(fij^)(3) = [(3 + v / 5)/2] N + 
[(3-V5)/2] N . 

In this space the translation operator shifting the system by two lattice sites can be diag- 
onalized yielding momentum eigenvalues P = I ■ 2tt/S where I = 0, ...,S — 1 (N = 2S). The 
energy eigenvalues of the Temperley-Lieb model in the sub-representation with fc-many singlets 
and momentum I ■ 2ir/S reference state are identical to those of the XX Z quantum chain with 
fc-many flipped spins (i.e. L — k up and k down spins with magnetization S = N/2) and twist 
angle 
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Note that the dependence of the eigenvalues of the XX Z reference system on ip is 7r-periodic. 
The quantum number I takes integer values for periodic and 7r-twisted boundary conditions of 
the Temperley-Lieb model, except for the 7r-twist and an odd number of fermions in the pseudo- 
vacuum. In such a case, the set of allowed values for I is shifted by 1/2 (Z = 1/2, 3/2, S — 1/2) 
which will not be pursued furthermore. 

The case S = (k = L/2, N = 0) is special. Here the relation between the twist angles is 
given by equating macroscopic loops 

1 T l + l = e ^ + e"^, (38) 

for periodic boundary conditions and 7r-twist. Explicitly we find 

ip = vr/3 resp. i In ((3 + y/E) /2) , (39) 

a real/imaginary value for periodic boundary conditions and 7r-twist. (For the general sl(n + l\n) 
case, the left hand side in (38) is to be replaced by 1 or d for periodic or 7r-twisted boundary 
conditions.) 

We read off the central charge from (28) with 7 = ir/3 and (39) as c = for periodic 
boundary conditions and c = 1 + (^ In [(3 + \/5)/2]) for the 7r-twist. The scaling dimensions 
corresponding to low-lying excitations in the sub-representation with fc-many singlets appended 
to a length N pseudo-vacuum (N = L — 2k) are obtained from (27) 

X = ^ l^(^) - £ i,o(^/3)] = ^ ([^(^) - ^,o(0)] - [£l(V3) - £i,o(0)]) (40) 

. . , 4S 2 + 9(m- tp/ir) 2 - 1 
= xs, m (<P) ~ *o,o(vr/3) = 12 ' (41) 

where the twist angle ip is as in (37) and the true ground-state corresponds to S = m = and 
twist angle in (39). 

Next we give the multiplicities of the momentum values in the space &P N of all pseudo- vacua 
of length N for periodic boundary conditions. First, the dimension of this space is (see appendix) 

-\ N 



d{N) :-.^ N) -.(i + ^ • (42, 

The subspace of Q P N of states with an odd number of fermions has dimension (see appendix) 

<f W = ^-cos(^). (43) 

The eigenstates of the translation operator by two sites in Q P N are constructed from the orbits 
with periods p (which are divisors of N) . Here we follow closely the treatment in [35] , the main 
difference is that we now deal with a translation operator by two lattice sites and hence an orbit 
with period p covers 2p sites. The dimensions of the orbits are denoted by a{p). Similary, we 
denote by cr (p) the dimension of the subspace of the orbit with period p containing states with 
an odd number of fermions. The dimensions of the orbits are related to the dimension of all 
pseudo-vacua by 

Y,°(p) = d(2S), (44) 
P\S 
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as well as 

^2 <y°{2 m p) = d°(2 m+1 M), (45) 

p\M 

where m and (odd) M are uniquely denned integers such that S = 2 m M. 
These equations are uniquely solved by 

'(SHE^ff) d ^ ( 46 ) 

r\S V ' 

and 

a°{2 m M) = E ^ (— ) d°{2 m+1 r) (47) 

r\M ^ r ' 

where [i is the Mobius function defined for integers q by 

1 , f or g = 1 , 
/•*(<?) := { {— l) s 5 if 9 is the product of s distinct primes, 
0, else. 

For periodic boundary conditions the multiplicity of the momentum in Q, P N is 
'2tt , \ - r 



M( ^j= e iKf (f + E §«•(?)■ («) 

7 r|gcd(S,0 V /J r|gcd(S,20 

reven, 2Z/rodd 



There are some interesting applications of the above results. The decay of 2-point functions 
at the critical point are governed by the lowest values of x with the right quantum numbers 
S, m, I, where, however, the application of selection rules is obscured as the space of pseudo- vacua 
il^r is rather high dimensional and 'violates' simple conservation numbers. 

With respect to the localization-delocalization exponent v the result 4/3 was found in [16] 
from the well established general RG scaling relation 

' (49) 



2-x' 



where the 'relevant' scaling dimension was identified as that of the 2-hull operator which corre- 
sponds to S = 2 and m = I = in (41). (For / = 1 the value of v is 1.7.) 

We like to note that the energy levels of the model by [15] show stronger logarithmic cor- 
rections, l/(L\nL), [17] than observed here, (lnL)/L 3 , cf. Sect. 3. This indicates that the two 
models correspond to rather different universality classes. 

Finally, we recall from [35] the finite temperature relation of a Temperley-Lieb vertex model 
with local Hilbert space of dimension d = 2n + 1 to the Heisenberg chain at same temperature 
plus an effective magnetic field h proportional to the temperature 

f(T) = fxxz(T,h), 2cosh^ = d. (50) 

In Fig. 5 we show results for the entropy and the specific heat of the first three cases of the 
Hamiltonians for n = 1,2,3. Note that the thermodynamical data at low-temperatures are very 
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similar for all cases with c(T) exhibiting a proportionality to temperature with slope ttc/3v with 
central charge c given by (28) and imaginary angle 92 with costp = d/2. At higher tempera- 
tures the curves deviate from each other. This is expected on physical grounds as in the high 
temperature regime the entropy has to approach Ind = ln(2n + 1). 




Figure 5: (a) Entropy S(T) and (b) specific heat c(T) data versus temperature T for the 
Temper ley-Lieb Hamiltonians with A = 1/2 corresponding to the sl(n + l\n) invariant net- 
work models with n = 1, 2, 3 (d = 3, 5, 7). 



If we impose 7r-twisted boundary conditions in imaginary time direction, i.e. if we define 
a pseudo-partition function by use of the super-trace, we find a vanishing free energy for all 
temperatures consistent with the low-temperature asymptotics with central charge c = 0. 



5 Conclusion 

In this paper we derived the complete low-lying energy spectrum of the quantum spin chain 
related to the sl(n + l|n) supersymmetric network of [16]. Our analysis consists of two major 
steps. 

First, we established a proper Temperley-Lieb equivalence of the quantum spin chain with 
periodic/twisted boundary conditions to a spin-1/2 XXZ chain with anisotropy A = 1/2 and 
suitable boundary conditions. We presented a representation theoretical approach as a simple 
alternative to the evaluation of the scaling dimensions and multiplicities in [33,34]. Second, the 
leading logarithmic corrections to the conformal spectrum were derived. Our results indicate 
that the models of [16] and [15] belong to different universality classes. 

The calculations are based on an asymptotic analysis of non-linear integral equations for 
the finite temperature and finite size properties of the Heisenberg spin chain. These equations 
are exact for arbitrary size. By these means, also the thermodynamical properties of the su- 
persymmetric Hamiltonians of Temperley-Lieb type were computed for arbitrary temperature. 
The specific heat exhibits a linear low-temperature behaviour with a slope corresponding to the 
central charge in the 7r-twisted sector. 
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Appendix 

In this appendix we give a short account of the derivation of the multiplicity formulas and a proof 
of the completeness of states. All formulae are to be understood as for the generic Temperley- 
Lieb representations in the limit of A = 2A — > 1 avoiding the explicit analysis of reducible but 
indecomposable representations. 

We begin with the characterization of the space of pseudo- vacua Q P N containing all states uj 
on a chain of length N with periodic boundary conditions annihilated by all operators 6i,...,6jv- 
If we label the canonical basis states of the local spaces V with j = 1, d (= 2n + 1) and define 
U : V — > V by Uij := Si^+i-j, the space Q 1 ^ is mapped by U <8> id <8> U <g) id <8> ... onto a space 
of states that are annihilated by the application of dual states (tp\ := Ylj=i 0>^+ 1 ~~ il an d 
(tjj := X^j=i( — 1) J+1 {j, d + 1 — j\ on (odd,even)- and (even,odd)-indexed nearest-neighbours, 
respectively. 

We introduce a ^-deformed version of this space by demanding annihilation by 
Ej=i q j ~ (d+1)/2 (j, d+l-j\ and £^ =1 (-l) i+ V W<m)/2 (j, d+l-j\. The dimension of the 
space is independent of q (see [35]) and can be easily enumerated in the limit (7—7-00. The 
condition for the states is simply that the sequence (d, 1) never occurs. The counting problem 
is now reduced to the evaluation of traces of powers of the d x d adjacency matrix A with 
entries 1 everywhere except for one at the position (d, 1). The only non-zero eigenvalues are 
a± = (d ± V d 2 — 4)/2, where a_ = \/a + . Hence 

d(N) := TtA n = a^ + a£, (51) 

which proves (36) for > 0. 

Defining the diagonal d x d matrix R := diag(+l, —1, +1, —1, ..) we find 

Tr(AR) N = d e (N) - d°(N), (52) 

where d e (N) and d°(N) are the dimensions of the subspaces with even and odd numbers of 
fermions, respectively. The only non-zero eigenvalues of AR are (1 ± i\/3)/2, hence d e (N) - 
d°(N) = 2 cos (f N) proving (43). 

Finally, we like to count the number of states covered by the constructed representations. 
The dimension of each /c-sector {k = 0, 1, L/2) is (^) and appears d(L — 2k) times for k < L/2 
and once for k = L/2. Hence the number of states is 

zVf) K 2fe+ ^ (L_2fc) ) + (^ 2 ) =j:{% L + - 2k =^+<) L =d L - (53) 

fc=o v y \ / / k=0 \ / 
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